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f^ , Abstract. The stochastic dynamics of heavy quarks in the fireball produced in heavy- ion 

I ' collisions is followed through numerical simulations based on the Langevin equation. The 

MhI modification of the final pr spectra [Raa) of c and b quarks, hadrons and single-electrons with 

respect to pp collisions is studied. The transport coefficients are evaluated treating separately 
the contribution of soft and hard collisions. The initial heavy-quark spectra are generated 
according to NLO-pQCD, accounting for nuclear effects through recent nPDFs. The evolution 
of the medium is obtained from the output of two hydro-codes (ideal and viscous). The heavy- 

^ ' quark fragmentation into hadrons and their final semileptonic decays are implemented according 

^^ , to up to date experimental data. A comparison with RHIC data for non-photonic electron 

spectra is given. 



o. 

^— ^ ■ 1. Introduction 

In a previous paper [Ij we evaluated the transport coefficients for heavy quarks in the QGP 
in the Hard Thermal Loop (HTL) approximation. The above coefficients were then used to 
^ ■ solve the relativistic Langevin equation in the case of a static medium, following the approach 

H I to thermal equilibrium. Here we improve our previous study in two important aspects. We 

provide a more accurate microscopic calculation of the transport coefficients and we perform the 
Langevin simulations for a fully realistic scenario, accounting for the hydrodynamical evolution 
of the medium. For independent studies based on the Langevin approach see [2l[3llll[5]. 

2. The Langevin equation in a dynamical medium 

The relativistic Langevin equation allows to study the stochastic dynamics of heavy quarks in 
the QGP. The algorithm to follow the evolution of momentum and position of the brownian 
particle is the following. We focus on a given quark at {xn,Pn) after n steps of evolution. We 
move to the local fluid rest-frame and update its position and momentum by the quantities 
Axn = iPn/Ep)Ai and 



^Pn = -VD{Pr.)fnAt + e{t)Ai^ -'nD{Pn)fnAi + g'\p^)rf {t)^ At, (1) 



where we take - in the fluid rest-frame - Ai = 0.02 fm/c. In the above we express the noise term 
through the tensor (we omit the "bar") 



g'' {p) ^ yMp)py + TMp) ('^''' - ff), (2) 

depending on the transverse/longitudinal momentum diffusion coefficients ht/l{p)^ ^ind the 
uncorrelated random variables rj^ , with {rf {t)r]^ {t')) = 5^^ du' ■ Hence, one simply needs to extract 
three random numbers rf from a gaussian distribution with 0" = 1. Concerning the problem of 
fixing the friction term r]Dip) and further details we refer the reader to Ref. [1]. We then go 
back to the Lab frame, getting the updated {xn+i,Pn+i)- The four- velocity and temperature 
fields u^{x) and T(x) of the background medium are obtained from the output of two different 
hydrodynamical codes [3 El EJ . 

3. Evaluation of the transport coefficients 

The coefficients Kt/l (p) reflect the transverse/longitudinal squared-momentum acquired through 
the collisions in the medium. Following [9] we introduce an intermediate cutoff |i|* ~ m'j-, 
{t = (P'—P)'^ = uj'^—q'^) separating hard and soft scatterings. The contribution of hard collisions 
{\t\ > \t\*) is evaluated in pQCD, from the diagrams Qq{q) — )■ Qq{q) and Qg — ;■ Qg and reads 
(employing the notation J^ = J dk/{2Tr)^ and Ptot = P+K—P' — K') 

and 

(g/q)hard. x 1 f ^B/f(^) /" 1 ± '^B/fI^') /" 1 ^,.,. |,|*wo x4r(4)/D ^ l^TJ |2 2 f a\ 

ip) = ^j^ ^T-A, — ^' — i ^^('*' - '*" )('") ^ ^^-*) l-^^/^l ^^' ^') 

where the squared amplitudes were evaluated in Ref. [lOJ- On the other hand in soft collisions 
{\t\ < \t\*) the exchanged gluon has a small virtuality and has thus "time" to feel the presence 
of the other particles. A resummation of medium effects is required and this is provided by the 
HTL approximation, leading to the compact formulas (x = uj/q) 
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and 



Cpg^ /-I*!* , r ^ |t|3/2 
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expressed in terms of the resummed gluon spectral function'p{\t\,x) = pL{\t\,x)+{v'^—x'^)pT{\t\,x) 
which can be found in [I] . Our results for k^/l (p) are shown in the left panel of Fig. [1] where a 
very mild dependence on the arbitrary intermediate cutoff is manifest. For large momenta kl 
largely exceeds kt- 

4. Numerical results 

For each explored case we generated an initial sample of 45 • 10^ cc and bb pairs, using the 
POWHEG code [11], with CTEQ6M PDFs. In the AA case we introduced nuclear effects in the 
PDFs according to the EPS09 scheme [12]; the quarks were then distributed in the transverse 
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Figure 1. Left panel: the momentum-diffusion coefficients kt/l{p) after summing the soft and 
hard contributions. The dependence on the intermediate cutoff |t|*~m|) is very mild. Right 
panel: the initial c and b spectra generated by POWHEG for pp and AA (with EPS09 nPDFs 
and nuclear /c^-broadening) collisions at -v/s^vAf ~ ^00 GeV. 
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Figure 2. Left panel: RAAiPx) for the c-quark spectra after the Langevin evolution in the 
QGP. The sensitivity on the hydrodynamical scenario (ideal vs viscous) appears to be small. 
On the other hand the dependence on the coupling g{fi) (at T = 200 MeV Os ~ 0.6/0.34 for 
fj, = 7tT/2ttT) is huge. Right panel: RAAiPr) for c, and b quark spectra, together with their 
weighted combination. 



plane according to the nuclear overlap function dN/dx±^TAB{x, y) = T^(x+6/2, y)TB{x—b/2, y) 
and given a further kx broadening on top of the "intrinsic" one. The initial pT spectra are 
displayed in the right panel of Fig. [TJ For each quark, after the initial free-streaming, at the 
proper-time r = t^ — z^ = ro we started following its Langevin dynamics until hadronization. The 
latter was modeled using a Peterson fragmentation function [13], setting e = 0.04 and 0.005 for c 
and b respectively; the hadron species were then assigned according to the branching fractions 
taken from |14U15j . Finally each hadron was forced to decay into electrons with PYTHIA [16j, 
using updated decay tables [T7j. The e-spectra from c and b were then combined with a weight 
given by the respective total production cross-section (c^c/66)- 

The results we present refer to minimum-bias Au-Au collisions at y/sj\f]y = 200 GeV. In Fig. [2] 
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Figure 3. Left panel: the effects of fragmentation and semileptonic decays, leading to a 
substantial quenching of the electron-i?/i^(pr). Right panel: Raa{pt) for the single-electron 
spectra from the decays of charm, bottom and their weighted combination. The latter can be 
compared with the experimental results by Phenix. 



we display the Raa{pt) = {dN / dtpx)^^ / {Nco\\){dN / dprY'^ for the final quark spectra. The 
dependence on the hydro scenario appears to be small. On the other hand the results display a 
huge sensitivity to the value of the coupling (hence, of the transport coefficients): a too strong 
coupling would overestimate the quenching effect. For the h spectra, the quenching induced 
by the medium is small: the most relevant effect is the larger production cross section in AA 
with respect to pp, arising from the anti-shadowing region in the nPDFs. In the left panel of 
Fig. [3] we consider the effects of fragmentation and decay into electrons: they both give rise to 
a quenching of the spectrum. Finally in the right panel we display, for a given scenario, the 
Raa{pt) of the single-electron spectra and we compare it to the results by Phenix |18j. Our 
results, obtained in a perturbative setup and referring to a quite moderate value of the coupling, 
appear in reasonable agreement with the experimental data. This also suggests that it should 
be worth addressing more carefully the role of collisional energy loss in the quenching of high-pj' 
spectra, so far mainly attributed to gluon radiation. 
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